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Abstract 

Let r be a nondegenerate geodesic in a compact Riemannian manifold 
M. We prove the existence of a partial foliation of a neighbourhood of F by 
CMC surfaces which are small perturbations of the geodesic tubes about 
r. There are gaps in this foliation, which correspond to a bifurcation 
phenomenon. Conversely, we also prove, under certain restrictions, that 
the existence of a partial CMC foliation of this type about a submanifold 
r of any dimension implies that F is minimal. 



1 Introduction 

Constant mean curvature hypersurfaces constitute a very important class of 
submanifolds in a compact Riemannian manifold {M"'^^,g). In this paper we 
are interested in families of such submanifolds, with mean curvature varying 
from one member of the family to another, which form (partial) foliations and 
which 'condense' to a submanifold F C A/ of codimension greater than 1. Our 
main results concern the existence of such families and, conversely, the geometric 
nature of the submanifolds F to which such families can condense. 

The simplest case, where F is a point, was considered by Ye a decade ago, 
|12|. He proved that if p E Af is a nondegenerate critical point of the 
scalar curvature function Rg , then there exists a neighborhood Li B p such that 
U \ {p} is foliated by constant mean curvature (for short CMC) spheres; in fact, 
the members of this family are small perturbations of the geodesic spheres of 
radius < p < po, and hence they have mean curvatures H — 1/ p ^ oo. 
Moreover, this foliation is essentially unique. Conversely, if a neighbourhood of 
p admits such a foliation, then necessarily Vi?g| = 0. In very closely related 
work, Ye ^HIj and by quite different methods (using inverse mean curvature 
flow) Huisken and Yau proved the existence of a unique foliation by CMC 
spheres near infinity in an asymptotically flat manifold (of nonnegative scalar 
curvature); this is of interest in general relativity. 

In this paper we study the existence of families of CMC hypersurfaces which 
converge to a (closed, embedded) submanifold F^ C M"+^, particularly in the 
case £ = 1. Define the geodesic tube 

Tp(F) {g e M"+i : dist<,((j, F) = p}; 
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this is a smooth hypcrsurfacc provided p is smaUcr than the radius of curvature 
of r, and we henceforth always tacitly assume that this is the case. The mean 
curvature of this tube satisfies 

n — f 

Hr,(T) = +0{l) as p \ 0, (1.1) 

np 

and hence it is plausible that we might be able to perturb this tube to a CMC 
hypersurface with H = {n — i)/{np). It turns out that this may not be possible 
for every (small) p > 0, or when i > 1, but we prove the : 

Theorem 1.1 Suppose that T is a simple closed embedded geodesic with nonde- 
generate Jacobi operator. Then there exist fco € N and sequences p'^. < p'l 0, 
for k > ko such that when p £ Ik := {p'k-.Pk), the geodesic tube 'Tp{T) may be 
perturbed to a CMC hypersurface Sp with H = . The Ep are nonintersecting 
and foliate the open set equal to their union, hence they form a partial foliation 
of some neighborhood ofT. 

The hypcrsurfacc Ep is a small perturbation of Tp{T) in the sense that it is 
the normal graph (for some function whose L°° norm is bounded by a constant 
times p^) over a small translate of T (by some translation whose L°° norm is 
bounded by a constant times p^), cf. §4 for the precise formulation. In addition, 
we have rather precise information on the location and width of the intervals 

pi-m^, = o{k-^i% ^^^^^ 

where A equals the length of T. Note that the nondegeneracy condition on F 
is a mild and generic one, and that there arc no stringent conditions on the 
curvature along F, as is the case in Ye's theorem when F is a point. 

The existence of gaps in this foliation stems from the fact that at certain 
radii, the Jacobi operators on the geodesic tubes 7^(F) become degenerate, 
and this substantially complicates certain analytic steps in the construction. 
However, this gap behaviour is a real phenomenon, and is linked to a bifurcation 
phenomenon, as we explain more carefully in §5. 

Recall that the index of a compact CMC or minimal submanifold is the 
number of negative eigenvalues of (the negative of) its Jacobi operator. We 
estimate the index of the leaves of the partial foliation we construct : 

Proposition 1.1 If p & Ik, then for k sufficiently large. 
Index (Sp) = Index (F) + 2 fc + 1. 



In particular, Index(Sp) — !■ oo as p ^ 0. It is easy to show that for a generic 
metric on M, the moduli space of CMC hypersurfaces condensing to F is a 
smooth one-dimensional set, i.e. a (possibly infinite) union of curves, and that 
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the index is constant on each component. Thus Proposition 11.11 shows that the 
leaves Ep for p £ Ik lie in different components of this moduli space. 

The final part of this paper concerns necessary conditions on T in order 
that a sequence of CMC surfaces condensing to F exists. We show that if there 
exists a sequence of CMC hypersurfaces S^-, each of which can be written (in 
an appropriate sense) as a normal graph over the geodesic tube Tp. (F) where 
pj 0, then F must be minimal. We defer to §5.2 for the precise statement of 
this result. We have proved this converse only under rather stringent conditions, 
but posit the following: 

Conjecture: Let T be a closed embedded l-dimensional submanifold of M and 
that there exist sequences p'^ < p'^ —^ and a partial foliation by CMC hyper- 
surfaces Sp, p€ h ■■= ip'k,Pk) ('W'i'th p'l < p'k^i), satisfying: 

(i) The mean curvature of Tip equals 

(ii) There exists a constant c > 0, independent of k and p € Ik such that 
Sp C {g e M"+i : distg{q,T) < cp}; 

(Hi) The norm of the second fundamental form of these hypersurfaces satisfy 

\Asp\ ^ c i for some constant c > 0, again independent of k and p Ik- 
Then T is a minimal submanifold. 

One might even be able to weak hypotheses (ii) and (iii) substantially, but 
even with these hypotheses, the proof is already probably difficult. We have 
chosen to prove this converse only under much stronger hypotheses in order 
to include one main calculation which explains why the minimality of F is the 
natural conclusion. 

On the other hand, our method of proof encounters serious analytic dif- 
ficulties when dimF > 1, and it is unclear whether there is a general result 
concerning existence of families of CMC hypersurfaces concentrating along a 
minimal submanifold of dimension greater than 1. The technical complications 
are due ultimately to the lack of sufficiently good estimates for (— Ar — A)^^ on 
Holder spaces when A ^ oo and A lies in a spectral gap. (When F is a curve, the 
spectral gaps are large and the Jacobi operator is an ODE, and such estimates 
are easy to obtain!). 

There are some parallels between our Theorem 1 and some recent results 
concerning solutions of the equation 

Au + f{u) = 

with vanishing Neumann data on a smooth, bounded domain VL C M". For ex- 
ample, the second author and Ritore ^ prove the existence of positive solutions 
to 

Au - + ii = 
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which concentrate along a minimal submanifold S as e \ 0. On the other hand, 
Malchiodi and Montenegro construct positive solutions of 

£^ Au + - w = 

which concentrate along dil as £ \ 0. As in the present paper, the same 
'spectral gap' phenomenon limits their results to domains f2 C M^. 

2 Geometry of tubes 

In this section we derive expansions as p \ for the metric, second fundamental 
form and mean curvature of the tubes Tp(r) and their perturbations. There is 
a famous and beautiful formula for the volume of these tubes, due originally to 
Herman Weyl, when AI has constant curvature, which has found applications in 
fields as diverse as geometric measure theory and statistics. We refer to Gray's 
monograph 131 for Weyl's formula and references to later work. 

2.1 Fermi coordinates and Taylor expansion of the metric 



We first consider the asymptotic development of the metric g in Fermi coor- 
dinates around F. This leads to an asymptotic formula for the metric on the 
geodesic tubes Tp(T). These computations are standard, and are described more 
systematically in 

Fix an arclength parametrization 'y{t) of F, t G [a, 6] :~ I, and denote by 
SNT the sphere bundle in A^F. Then 



is a diffeomorphism when p is small enough. Choose a parallel orthonormal 
frame £'i,-- - ,£"„ for A^F (along (a, 5), say). This determines a coordinate 
system 

X := {xo,Xi ...,Xn) I — > exp^(^^)(a;ii;i + . . . + a;„£'„) F{x), 

and the corresponding coordinate vector fields Xa ■= F^,{dx^)- We write x' = 
{xi, . . . ,Xn), and adopt the convention that indices i,j,k,... € {!,..., n}, 
whereas a, /3, ... € {0, ... , n}. 

Remark 2.1 For simplicity, we identify the metric g on M and its pullhack F*g 
on some neighbourhood m M x R". This allows us to use the linear operations 
on the latter space. With slight abuse of notations, we identify F{xo,x') with 
(xq, x') and Xa with dx^ ■ 

We also use cylindrical Fermi coordinates. Thus let r = xf + . . . + x"^, 
which by Gauss' lemma is the geodesic distance from x to F. The vector 



near F 



SNT 3 {t,v) ^ exp^(,)(H e Tp{T) 




(2.3) 
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is the unit normal to the geodesic tubes. 

We have arranged that the metric coefficients gap — {Xa,Xp) equal Saf3 
along r. We now compute higher terms in the Taylor expansions of these func- 
tions. In the following, the notation ©(r™) indicates a function / such that 
it and its partial derivatives of any order, with respect to the vector fields Xq 
and Xi Xj, are bounded by Cr"* in some fixed Tpg{r). Also, we shall compute 
the metric coefficients at a point q := F{xo,x') in terms of geometric data at 
p := F{xo,0) and the radius r = q). 

We begin with the expansion of the covariant derivative: 

Lemma 2.1 For a, f3 = 0, . . . ,n, 

n 

V;f„X^ = ^0(r)X^, (2.4) 

7=0 

and for a ~ f3 ^ 0, we record the more precise expansion 

n n 

Vxo Xo^-J2 W^j'^o) ^o)p X, X, + O^'^'') ^7- (2-5) 



Proof: Anywhere on F, 

Vxo^o = ^XoX, - Vx.Xo = Vx^X, = 0. 

The vanishing of the first two terms is obvious since F is a geodesic and the Xi 
are parallel along it. Because we are using coordinate vector fields, Vx^Xfj = 
^XfjXa for any a,/?, even away from F, and this implies the vanishing of the 
third term. Since any X g NpT is tangent to the geodesic expp(sAr), and so 
^Xi+Xj {Xi + Xj) = at p, hence \7xiXj + \/ XjXi — there. Combined with 
the symmetry statement, we obtain that the final term also vanishes. This now 
gives (|2.4|l . 

Next, using H2.4|l we get 

X,{VxoXo, X,) = {Vx^XoXo, X,) + {VxoXo, Vx^X,) 
= {Vx,Vx,Xo,X,)+0{r-') 

= {R{X,,Xo)Xo,Xj)p + {Vx„yx,Xo,X,)p + 0{r^) 
= {R{X,,Xo)Xo,Xj)p + 0{r) 

This implies 1^3)) . □ 

Our next result gives the expansion of the metric coefficients in Fermi coor- 
dinates. The expansion of the gtj, i,j = 1, . . . , n, agrees with the well known 
expansion for the metric in normal coordinates, cf. |10|. [B] or jl5| . but we 
briefly recall the proof here for completeness. 
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Proposition 2.1 In the same notation as above, we have 

g^M - + i {R{Xk,Xi) Xi, Xj)p Xk xt + 0(r3) 

5o.(g) = 0{r^) (2.6) 

5oo(<z) = l + (i?(Xfc,Xo)X,,Xo)pXfea;, + 0(r3). 

Proof: The function 

Xkgap = ^ X^Xa-iXp) + {Xa,y X^Xji) 

vanishes on F, and thus the first order terms vanish in all of these Taylor ex- 
pansions. 

To compute the second order terms, it suffices to compute 

Xk Xk Qap {p) = Xk Xk {Xa , Xp) 

= {\/l^X^,X0) + {X^,\/l^Xi,)+2{\/x,X^,\/x,X0) 

and then polarize (i.e. replace Xk by Xk + Xg, etc.). By (12.4(1 . the final term 
vanishes. Also, 

Vj,^X^=Vx,yx^Xk^Vx^Vx,Xk + R{Xk,X„)Xk. 

First let a = J > f and compute the first term on the right. Since 
^ x^ xX = on F for any X which is a constant linear combination of the 
Xi, we have 

= '^Xk+eXj'^Xk+eXj (Xk + sXj); 
equating the coefficient of e to gives V x,'^ XkXk = —2Vxk'^Xk^j: ^^^d hence 

'^^'xk'^O — R{^k,Xj) Xk, 

so finally 

2 

Xk Xk gij = - {R{Xk,Xi) Xk-Xj). 

The formula for the second order Taylor coefficient for gij now follows at once. 
When a = 0, Vxo^Xj^Xk = on F, so 

^Xk-^o = R{Xk,Xf)) Xk- 

from which it follows that 

Xk Xk goo = 2 {R{Xk, Xq) Xk,Xo) 

and this gives the formula for 500 • 

The second order Taylor coefficient for goi has not been given because it is 
not needed later. □ 



6 



2.2 Perturbed tubes and their mean curvature 

Wc now describe a suitable class of deformations of the geodesic tubes Tp{r), 
depending on a section $ of NT and a scalar function w on the spherical normal 
bundle SNT. One of the main technical parts of this paper, which occupies the 
rest of §2, is the computation of the mean curvature of these hypersurfaces, at 
least asymptotically in p and for sufficiently small $ and w. 
The spherical normal bundle is locally trivialized by the map 

(a, 6) X ^"-1 9 {xo,e) ^ {j{xo),J^0jEj) e SNT. 

Fix p > 0, and define 

G{xo,e) := F {xo, p{l + w{xo, 6)) 6 + $(a;o)) ; 

the image of this map will be called Tp{w, Thus Tp{w, $) is obtained by first 
taking the normal graph of the function p w over the tube of radius p in NT 
and then translating by $. In particular 

rp(o,o) = rp(r). 

It will sometimes be useful to calculate using a coordinate system 

with associated coordinates vector fields Yj = dyjT. In particular, we regard G 
as a function of {xq, y) and write 

G{xo,y) := F{xo,p{l + w{xo,y)) T + $(a;o)) . 

Two different types of Holder spaces will be used to measure regularity 
of fimctions on SNT and sections of NT: first, we use the ordinary Holder 
spaces C'"'"(5'7Vr), C™'"(r,7Vr), but we shall also use modified Holder spaces 
C^'°'{SNT), C™'"(r,iVr) which are based on differentiations with respect to 
the vector fields pdx„ and dy^ (where y is any local coordinate system on S"""^, 
see above) . Note that this is tantamount to using the rescaled variable s = xo/p 
since ds = p 9xo • We shall assume that 

n 

$(xo) = Y,4>MXj e c2'"(r,ivr), w e cI'''{snt). 

i=i 

For p S r, let Sp denote the spherical fibre of SNT over p. Any function w 
on SNT decomposes into a sum of three terms 

W = Wq -\- W + W, 

where the restriction to any Sp of each of these terms lies in the span of the 
cigcnfunctions ipj{d) on S'^~^ with j = 0, j = 1, . . . , n, and j > n, respectively. 
The first component, wq, is a function on F itself. Next, the eigenfunctions 
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fji 1 j ^ ^''c the restrictions to 5"^^ of linear functions on M", so 
any linear combination of them can be identified with a translation in R" (the 
linear function x ^ a ■ x being identified with the translation x x + a). 
Correspondingly, the summand it) is canonically associated to a section $ of the 
normal bundle NT. 

We shall typically assume that the functions w has 'linear component' it) = 0, 
and shall regard the linear part of the perturbation as a section of AT, as just 
described. 

It will be fundamental in the analysis below to regard w as a, function of 

s = xq/p and yj, but $ and 7 as functions of xq (in particular, whenever we 
write we mean dxQ^ ^OxQipji^o) Xj)- However, we sometimes also write 
G = G{s,y). For example, the tangent space to 7^(w,$) is spanned by the 
vector fields 

Zo = G,{d,) = p{Xo + d,wT + ^'), 

(2.7) 

Zj = G^dy^) = p{{l + w)Yj +dy^wT), j = l,...,n. 

Definition 2.1 In the following, L{w, $) denotes any expression which is a 
linear differential operator ( of order at most 2) in w and $ which satisfies 

\\L{w, $)||co.» < c (||w||c^,-(s^r) + Il*llc2-(r,jvr)) , (2.8) 

where c is independent of p. Similarly, Q{w, $) denotes any nonlinear differ- 
ential operator ( of order less than or equal to 2) in w and $ which vanishes 
quadratically in the pair {w, $) and such that 

||Q(w2,$2) - g(wi,$i)||co,. < c ^sup^ (\\wi\\cf-(SNr) + ll*»llc=.-(r,wr)) 
X ^||«;2 - «;i||c2,a(5;vr) + 11*2 - $i||c2><»(r,ivr)^ (2.9) 

Here the spaces C°'" are either equal to C°'"(S'iVr) or C°'"(r, iVF) according to 
the range of L and Q. Finally, terms denoted 0{p^) are hounded inC"^'"'{SNT) 
or ^"'"(r, AT) hy C , where the constant C does not depend on {w, $) or p. 

2.3 The first fundamental form 

The next step is the computation of the coefficients of the first fundamental 
form of Tp{'w,^) with respect to the coordinates {s,y). At the point 

q = F{ps, p{l + w{s, y))T{y) + $(ps)) 
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(and p = F{ps, 0)), we obtain directly from 1)2. 6|l that 

{Xo,Xo)g = l + 0(p2)+pL(u.,$) + g(«;,$) 
{X,,X,)g = (X„X,>p + 4(i?(T,X,)T,X,)p + 0(p3) 

+ f [{R{T, X,) X,)p + X,) T, (2.10) 

+ p2L(w,$) + Q(w,$) 

(X„Xo), = 0(p2)+pL(u;,$)+Q(w;,$). 

We use these expansions to obtain the expansion of the first fundamental 
form of Tp{^, w). 

Proposition 2.2 We have 

p-^Zn,Zo) = l + Oip^) + pL{w,<i>)+Q{w,^) 

p-^Zo,Zj) = 0(p2) + L(w,$) + Q(u;,$) 

p-HZ.,Z,) = {Y,,Y,)p + i{R{^,Y.,)T,Y,), + 0{p^) 

+ 2 w (r„ y,)p + f [(i?(T, y,) y,)^ + {r{t, y^) $, 

p2 L{w. $) + (5(u;, $). 

(2.11) 

Proof : The first equation is clear. We give more details about how to derive 
the second and third estimates since the same argument will be used frequently. 
First, it follows from (|2.1()|l that 

(T,r,), = (T,Y,)p + 4(i?(T,T)T,y,)p + 0(p3) 

+ f [(i?(T, T) + (i?($, T) T, 

+ p^L{'w,'i>) + Q(w,$) 

However, when w = $ = 0, (T,Yj)g = since T is normal and Yj is tangent 
to 7^(0, 0) then, so that the sum of the first three terms on the right, which is 
independent of w and must also vanish. This, together with the fact that 
i?(T, T) = implies that 

(T, Y,), = I (i?($, T) T, + p2 ^(^^ ^) ^ (2.12) 

In particular, we get {T,Yj)q — pL{w,^) + Q(w,$)). The second and third 
equations follow directly from this. □ 
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2.4 The normal vector field 

The next task is to find expansions for the unit normal to Tp{w, <J>). We begin 
with the preparatory 

Lemma 2.2 The following expansions hold 
(T,T), - l + p^L{w,<i>)+Q{w,^) 

{T,Zo)q = pLiw,^)+pQiw,<l>) 

{T,Zj), - pdy^w + !^{R{<S>,T)T,Yj)p + p^L{w,<^)+pQ{w,<i>) 

Proof : These follow from (|2.1U|) . As at the end of the last subsection, we are 
using that (T, T)g = 1 when w — ^ — Q and i?(T, T) = to obtain the first two 
expansion. The second expansion follows from the fact that (T, Zo)^ = when 
w = $ = 0. Finally, to obtain the last expansion, we use that {T,Zj)q — 
when w = <I> = as well as the first expansion and H2.12|l . □ 
We can now proceed with the expansion of the unit normal vector field to 
Tp(^«,$). 

Proposition 2.3 The normal vector field N to 7^(w, has the expansion 
N -T + j:';ila,Y, + {L{w,^)+Q{w,<i>))Xo 

(2.13) 

where the coefficients aj are solutions of the system 

n-1 

J2 «j , Yi}p = dy^w + ^(i?(<f>, T) T, Y^p. 



Proof : Define the vector field 

n-1 

iV := - T + ao ^0 + ^ , 

by choosing the coefficients aa so that that N is orthogonal to all of the Z^- 

It follows at once from Lemma and (|2.11() that pua = L{w,^) + Q{w,^) 
for every a. Plugging this back into each of the equations {N,Zi)q = (thus 
neglecting the orthogonality condition when a = now), we find that the aj 
are solutions of the system 

n-1 

{Yj,Y,)p ^-dy^w+ -(i?($, T) T, Yi}p + pL{w, $) + - Q{w, $). 

Recall also that Zj ~ pYj + pL{w,^) so that OiZi — paiYi + pL{w,^). Finally, 
we have 

\N\g = l + p2L(u;,$)+Q(u;,$). 
This gives (ITT^ . □ 
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2.5 The second fundamental form 

The most arduous step is the computation of the second fundamental form. 
To simphfy the computations below, we henceforth assume that, at the point 
T{y) G S^-\ 

(r„ Y,)p = (5y and Vy.Fj = 0, i, j = 1, . . . , n - 1 (2.14) 
(where V is the connection on TS"^^). 
Proposition 2.4 The following expansions hold 

p-^N,Vz,Zo)q = p{R{r,Xo)T,Xo)p + Oip^) 

- i a> - ($", T>p + (i?,(T, Xo) ^, Xo)p (2.15) 
pL(«;,$) + ig(«;,$), 

p-^ {N, VzoZj)q - 0{p) + i L{w, $) + i Q{w, $), (2.16) 

p-^N,^z,z,), = i + |p(i?(T,y,)T,r,)p + o(p2) 

+ pL(u;,$) + iQ(u;,$) 
p-2(iV,Vz.^,), = 0(p) + ii(u;,<i>) + ig(«;,$), i ^ j. (2.18) 

Proof : First note that by Lemma ITTl 

n 

Vx„ Xp\^ = {0{p) + L{w, $) + Q{w, $)) X^, (2.19) 

7=0 

since the coordinates of q depend on w and $. Hence , as 



Vz^Xp = Y,{0{p^)+pL{w,^)+pQ[w,^)) X^, (2.20) 

7=0 

which follows from H2.19|l and the fact that Za= p ^^{'^ + L{'w, <&)) X^. 
We will also use that 

n 

N + T = Y,{L{w,<^>) + Q{w,<^)) X^, (2.21) 

which follows from H2.13|l . Finally, we will need the expansions 
(T,Xo), = pL(ii;,$)+Q(u;,$) 
(T,r,), = pL{w,<^)+Q{w,^) 

11 



(2.22) 



whose proof can be obtained, as in Lemma [2. 21 starting from (|2.1U|) . 
Estimate (|2.15|) : We must expand 

(TV, Vzo^o), = P'^ ((iV, Vzo^o), + {N.VzAdswT)), + (TV, Vz„$%) 

The estimate is broken into three steps: 
Step 1 From H2.13|l and Lemma [2. 21 we get 

(7V,T), = -(T,T>, + E;=i«,(iG-,T), + (L(u;,$)+g(u;,$))(Xo,T), 

+ (p2LK$)+Q(z„,<i>))(X„T), 

= -1 + L{w,^) + Q{w,^) 

Substituting N ^ -T + N + T gives 

{N, Vz„T), = a,(T, T), + (TV + T, Vz,T),; 

by Lemma 

93(T,T), =p2^(u;,$)+Q(u;,$), 
and and itT^ imply 

(iV + T, VzoT), - p2 ^(w;, $) + p Q(«;, 
CoUecting these estimates we get 

(iV,Vz„T), =p2L(^,$) + g(^,$). 

Hence we conclude that 

{N,Vzo{dswr))q = dlw{N,T)^ + dsw{N,Vzo^)i 
= -d^w + Q{w,^) 

Step 2 Next, 

n 

(TV, Vz„$'), = P (A^, + E ^-o<^^- ^.>9 
where $"(a;o) := c^xo'?^i(^o) ^j - From H2.20|l . we have 

Also, using the same decomposition of N, and employing H2.21|l and H2.10|l . 

(TV,$">, = -(T,$"), + (TV + T,$'% 

= -{T,<i>")p + p^L{w,<f) + Q{w,<^) 
12 



Collecting these gives 

(iV, Vzo*% - (T, $">p + Liw, $) + pQ(w;, $). 
Step 3 Expanding Zq gives 

n 

(TV, Vzo^o), = P {N, VxoXo)g + pdsw {N, VrXo)g + pY.^o {N, Vx.^o), 

(2.23) 

With the help of H2.19|l and (|2.21|) , we evaluate 

(A^,VtXo>9 - 0{p)+L{w,<^) + Q{w,<i>) 
{N,Vx,Xo), - 0(p)+L(w,$) + g(?«,$) 
(A^ + T,Vxo^o), = pL(«i,$) + Q(u;,$), 
and plugging these into (|2.23|) already gives 

{N, VzoXo)g = -p (T, Vxo^o)g + p' i(^«, + pQ(w, $) 
Using (|2.5|) in Lemma [2. II we get the expansion 

Vxo^ol, = -T.%iP{R{X„X^)T,X^)^pX,+0{p^) 

- Y.%i{R{Xj.Xo) $, Xo)p X, + pL(7«, $) + Q{w, 

and so 

(7V,Vz„Xo), = p2 ^;^^(^(^^.^^^)Y^x^)^(T,X,),+0(p3) 
+ P E;'=iW^..^o)$,^o)p (T,X,), 

Finally, using H2.10|l again, we conclude that 

(iV,VzoXo), = p2(i?(T,Xo)T,Xo)p + 0(p3)+p(i?(T,Xo)$,Xo)p 

+ /72l(u),<i>) + pg(w,<i>), 

which, together with the results of Step 1 and Step 2, completes the proof of 

the first estimate. 

Estimate (|2.16|) : Decompose 

{N,Vz„Zj)q = p{N,Y,)qdsW + p{N,T)qdsdy^w 

+ p{l + w) (iV,Vz„rj>g + p(A^,VzoT),. 

As above we use (|2.21|) and (I2.22|l to estimate 

(iV, Y,), = -(T, + (TV + T, = L{w, $) + Q{w, $) 

13 



Similarly, by Lemma IT^ and H2.21|l . 

(iV,T), = -l + L(^i;,$)+Q(^i;,$) 
But now, by H2.21|l and (|2.2U|) . we have 

(TV, Vzo r,), = O(p') + + pQ{w, $) 

and 

(TV, Vzo T), - + $) + p Q(«;, $), 

and the proof of the estimate follows directly. 

Estimates (|2.17|) and (|2.18|) : Observe that, thanks to the result of Proposi- 
tion we can also write 

1 " 

where 

n 

N = {L{w, $) + Q{w, $)) Xq + J2{P^ L{w, $) + Q{w, $)) Xj. (2.24) 
Now write 

{N,\/z,Z,,), = (iV,Vz^.,Z,% 

= \ ((Vz,T,Z,.), + (Vz^.,T,Z,%) 

- T-p ELi ((Vi?, (a. ^0, Zy), + (Vz,, (a. Z,), Z,)^ 
+ i ((iV,Vz^Z,.), + (iV,Vz^,Z,%) 
-i (9,^. (TV,Z,v)|,+9,^., (iV,Z,)|,) 

Step 1 By (|2.19|l . we can estimate 

Vz, ^j' P dy, wYj> + pdy^dy^,wT 

+ pil + w)Vz, Yr + P dy^, w Vz, T 

+ ELi W) + pLiw, •f) + p2 Q{w, $)) Xfe, 

Observe that the coefficient of Xq is slightly better than the coefficient of the 
other Xk since the first two terms only involve the Xk ■ Using this together with 
H2.24|l we conclude that 

(TV, Vz,Z,,)g + (TV, Vz^,Z,)g = p3 ^) +pQiw, $) 
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step 2 Next, using (|2.24|) together with H2.10|l . we find that 

dy^ {N, Zj,)\g + dy^, {N, Z,> I, = p3 ^) ^ ^ Q(^^ ^) 

Step 3 We now estimate 

It is convenient to define 

YT^ ((Vz, (1 + ^i;) T, Z,,)g + {\/z, (1 + w) T, Z,),) , 
It follows fi'om Lemma [2.21 that 

hence it is enough to focuss on the estimate of A'^^,. To analyze this term, let 
us revert for the moment and regard w and $ as functions of the coordinates 
{t^y) (rather than (s,y)), and also consider p as a variable instead of just a 
parameter. Thus we consider 

F{p, t, y) = F{t, p(l + w{t, y))T{y) + <i>{t)) . 

The coordinate vector fields Zj are still equal to {dy- ) , but now we also have 
(1 + w)T — F^,{dp), which is the identity we wish to use below. Now, from 
H2.11|l . we write 

4/ = {{^o,Z,,Z,,), + {Va^Z,,,Z,),) = ^ dp{Z,,Zy)\, 

Therefore, it follows from (|2.11|) in Proposition 12. 21 that 

= T^dp[p^Y,,Y,,), + ^{R{T,Y,)T,Y,,), + 0{p^) 

+ 2p^w {Y,,Y,,)p + ^ ((i?(T, Y,) Y,,)p + (i?(T, Y,,) $, Y,)^) 
+ p'^ L{w,^) + p^ Q{w,^)] + pQ{w,<i>) 
= ^ [2p{Y,,Y,,)p + I p3 (^(T^ Y,) T, + 

ipw (r„ ry)p + p' ((i?(T, Y,) Y,,)p + (i?(T, Y,-o 

+ +pg(w,<I>)] 

+ 2pw{Y„Y,,)^ + p2 ((i?(T, y,) $, y,,)p + (^(T, r,') $, y,),) 

+ p^ L{w,^) + pQ{w,^) 
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step 4 Finally, we must compute 

Bjj' ■■= {^z,{a^ Zi), Zj,)q + (Vz^, {a^ Zi), Z^), 
= (Zi, Zj,)q dy^ai + {Z,, Zj)qdy.,ai 

+ a, {{Vz,Z,, Zj), + (Vz^,Z„ Z,),) 
= {Zi, Zji)q dy^ai + {Zi, Zj)q dy^,ai 
+ a, {{\7z,Z,, Zj), + (Vz.Zy, 

= {Zi, Zji)q dy^ai + {Zi, Zj)q dy^,ai + Oli % {Z^,Zy)q 

Observe that, H2.14|l implies 

^y^{Y^,Y,.)y = ^. 

Using this together with (|2.11|) and the expression for the given in Proposi- 
tion we get 

a,dy^ {Zj,Z\)q^ L{w,^) + Q{w,^) 
It follows from (|2.11() and the definition of again that 

{Z,, Zy)qdy^a, = p2 (K„ y^.,)p dy^a, + L{w, $) + p^ Q{w, $) 
Therefore, it remains to estimate {Yi,Yji)pdy-ai. By definition, we have 



J2 a. {Y,Yy)p = dy^,w + f (i?(<i>, T) T, Y,,)p 



Ui [U, I]'jp = Uy.,W -h ^ 

i=l 

Differentiating with respect to yj we get 



E ((^'' dy.^^ + {Y^,Y,,)p) = dy^dy^,W+ ^ dy^ (i?(<i>, T) T, Y,,)j> 



P 

3 

(2.25) 

Again, it follows from H2.14|l that dy.(Yi,Yji)p = 0. Moreover, using H2.2()|l . we 
first estimate 

Vz, T = + 0{/) + pL{w, $) + pQ(w;, $); 
and, using in addition H2.I4|I . we also get 

Vz,r/ - a T + 0(p2) + pL(u;, $) + pg(w, $) 
for some a G M. Reinserting this in (|2.25() yields 

Eti {Y^ - Yr)p dy, = dy^ dy^, w + f ) T, y,, ), 

+ f T) y, , yy )p + p' L{w, $) + p2 Q(y;, $), 
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since i?(T, T) = 0. 

Collecting these estimates, we conclude that 

Bj, = d^w + ^ {R{<P, Y,) T, Y,)p + p4 L{w, $) + Q{w, $) 

since T) Yj,Yj)p = and also that 

Bjj, = p2 L{w, $) + p2 Q{w, 

when j ^ f . With the estimates of the previous steps, this finishes the proof of 
the last two estimates! □ 

2.6 The mean curvature 

Collecting the estimates of the last subsection and taking the trace, we have now 
proved that the mean curvature H{w, $) of the hypersurface Tp{'w, $) satisfies 

npH{w, $) - (n - 1) = 
(f (i?(T, Xo) T, Xo)p - i Ric(T, T) ) p2 + o{p^) 
- (a> + As,.-iw + {n-l)w)-p ($" + Xo) Xo, T)p (2.26) 
+p'^L{w,^)+Q{'w,^). 

(We recall that if is an orthonormal basis of TpM, then 

n 

Ric(T, T) - ^ (i?(T, ^„)T, 

a=0 

Denote by (Xj^ipj) the eigendata of A^n-i, where the eigenfunctions are 
orthonormal and counted with multiplicity. 

A most important observation is that the second and third terms in the 
expansion oinpH are quadratic in the coordinates Xj. Hence, when — w — 
we have 

{npH - {n-l),Lpj)^. = 0{p^), j = l,...,n, 

(2.27) 

(npi/- (n- 1),(^j)l2(S"-i) = C'(p2), j7^1,...,n, 

or in other words, writing / — npH — (n — 1), then / = 0{p'^) but its ^^(5*""^) 
projection over i^si, . . . , satisfies f = 0{p^). 

3 Jacobi operators 

In this section we examine the mapping properties of some linear operators 
which appear in the expression of the mean curvature of 7^(u', $) given in (|2.26() . 
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3.1 Definitions 

The two linear operators appearing in the third hnc of 12.26|l are 
w I — > Csnfw := d'^ w + Ae w + {n — 1) w, 

(3.28) 

The latter is the Jacobi operator on T corresponding to the second variation of 
the length functional on curves, while (up to a multiplicative factor) the former 
is the Jacobi operator for the second variation of the area functional about a 
Euclidean cylinder M x S'"'"^(p). 

Recall that the geodesic F is said to be nondegenerate when J is invertible, 
i.e. if the equation 3 $ = has no nontrivial solutions on all of F. For a generic 
metric on M, it is well known that all closed geodesies are nondegenerate. 

On the other hand, since it is already naturally expressed in terms of the 
scaled coordinate s, 

CsNF ■■ C^^'^iSNT) C°'"(5iVF) 

is bounded uniformly in p. We can analyze this operator using the eigendecom- 
position for Ag on S"^^. As in §2.2, if the eigenfunction decomposition of w is 
given by 

w{s,e) = w-iis) ip-iie), 

j>0 

then w decomposes as wq + w + w, where 

n 

w := Wj (fij and w Wj ipj . 

j=l j>n 

We denote by XIq, 11 and H the projections on to these three components, 
respectively. From now on, we assume that we are working with functions w 
such that n It; = 0, and thus we only need to be concerned with the operators 
{CsNr)o and Csnt induced on the two other components. Note in particular 
that 

3.2 Mapping properties 

We now study the mapping properties of -3 and (the components of) Csnf- 
We first note that 

3 : C2'"(F, NT) — > C°'"(F, NT) 

is an isomorphism when F is a nondegenerate geodesic. 
Next, we also assert that 

CsNT ■■ nc2'"(5iVF) — > UC^p^iSNT) 
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is an isomorphism with inverse uniformly bounded as p ^ 0; this follows from 
the fact that A^i + (n — 1) < — C < on this subspace. Details are left to the 
reader. 

Finally, it is clear that 

(£sJvr)o : noC^'«(5iVr) noC^p'"{SNr) 

is bounded for every p > 0, but is only invertible when 

Vn - 1- 4 27rZ; 
P 

in the exceptional cases, there is a two-dimensional nullspace spanned by 



cos(vn^^s), sin(vn — Is), 
and hence a two-dimensional cokernel. To determine the norm of its inverse 



when \/n — 1 A/p ^ 2 ttZ, suppose that {£sNr)o v = f. Then 

\/n — l v{s) = sin(\/ri^~Ts) (a + /g cos(\/« — 1 o") f{a) da) 
— cos{y/n — Is) (/3 + Jq sin(v'n — la) f{a) da) 
where the constants a, /3 are chosen so that w is A/p-periodic. We find 

p (1 - cos(vn - lA/p)) 
for some constant c > 0, independent of p, and from this we have 

(3.29) 

again for some constant c > independent of p. Note that when / S C^, there 
is an equivalent formula 

{n — 1) v{s) — f {s) = s,m.{\/n — Is) (a — sm.{\/n — la) dsf{a) da) 

+ cos{\/n — Is) (/? — cos{\/n —la) dsf{a) da) , 

where again a, (3 are chosen so that w is A/p - periodic. We now obtain 

p (1 — cos(v?^ — lA/p)) 
for some constant c > independent of p, so that 



rllcp'^csivr) - 



( \\^\h-iSNT) + p(i_cos(Vn-lA/p)) l'^-^"^°°("^")0 ■ 

(3.30) 
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4 The constant mean curvature foliation 



We now use the results of §2 and §3 to perturb Tp{r) to a constant mean 
curvature hypersurface, at least for p sufficiently far from values where {CsNr)o 
is degenerate. 

According to the analysis of §2, we must find w E Cp'°'{SNr) and $ G 
C2>"(r,iVr) such that 

= n- 1 (4.31) 

Let us denote by 

/ := ^ {R{T, Xo)T, Xo)p ~ ^ Ric(T, T) + 0{p^), 

the inhomogeneous term appearing in H2.26(l which corresponds to the mean 
curvature when w — ^ — 0. As usual, this decomposes into three components, 
/o + / + /, where / corresponds to a section of the normal bundle which we 
write as p^. We are searching for w = wq + w and <& which satisfy the coupled 
system 

{CsNr)o Wo = fo + P^ L{w, $) + Q{w, $) 

a$ = * + + (4.32) 

CsNrw - f + p^L{w,<P) + Q{w,^) 
We use the function space 

:= Uq Cl^°'{SNr) e C2'"(r, NT) QllCf^iSNT), 

where, for S — {wq, $,i&), 

||S||£2.„ := (l-cos(Vn- lA/p)) ||wo|lc2-°(SArr) + ll'^llc2,=(r,Arr) + \\w\\cl-'{SNr)- 

The linear operators appearing on the left in (|4.32l) are all invertible provided 
y^n — 1 A/p ^ Z. Thus, multiplying by their inverses, we rewrite this system as 

and so we solve our problem by finding a fixed point of 91 in £p'". 
Lemma 4.1 Write Eq := *Tt(0). Then for \Jn — ^ ^/ p ^ we have 

l|2o||,..<f P^ 

for some cq > 0. 
Proof : Clearly 

ll/ollc;^-°(SAfr) + P'^ WdsfoWh-^iSNV) < cp^; 
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moreover 

ll/llc°'°(swr) ^ and ||^'||cO'°(r,A'r) < cp^. 

Since, by definition, 

So = {{CsNr)o^ fo,d^^ *o, [CsNry^ /), 

the result follows from (|3.30() and the uniform boundedness of the inverses of 
the inverses of these linear operators as p ^ 0. □ 

Next, from the properties of the operators L and Q we deduce the 

Lemma 4.2 There exists a constant c > such that, for the same Cq as in the 
previous Lemma, and for any Si,S2 G £p'°' satisfying 

||Si||g2,o < cop^ 



we have 



||m(S2)-m(Si)L.,„ <c- 



^"^ ~ (1 - cos(Vn^A/p))2 
Proof : It follows from H2.8|l that 



F2 - =1 k^.' 



" 1 — cos(v?^ — lA/p) 



.2 



if w = wq + w and ^ — (wq, w). Moreover, if ||Si||g2,Q < co p^, then we have 
from lEini) 

\\Q{W2, $2) - QK, < C- ^L=— — IIS2 - Sill 2,„. 

(1 — cos(v7i — lA/p))^ '' 

where Wi = wo,i + Wi and S = (wo.i, $i, Wi). Now, the result follows at once 
from the inequality 

||*||c"-°(r.Arr) < cp"" ||$||co.»(r,Arr) 

and the uniform bounds on Cgl^Y" ^^'^ 5"^, and the bound on {CsNr)o^ given 
in (|3.29|l . Details are left to the reader. □ 

Collecting these results, we now have the 

Proposition 4.1 Fix a £ (0,1). Then there exists a Ci > such that if k is 
sufficiently large and p satisfies 

\ Cl 27: 1 Cl 

^ ; T . P ^ 



k + 1 /c9/4 - V^T^ A - k fc9/4 

then there exists a solution {wq,w,^) of {4--^'/^ in fp'". This solution satisfies 

\\{W0,W,^)\\^2.a^ < COP^ 
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Proof : It is easy to check that 



(1 - cos(Vn^A//9))2 

is as smah as we want, provided ci is chosen large enough. It is then easy to 
check that, when k is large enough, *yi is a contraction from the ball of radius 
Co into itself. □ 

This proposition yields the existence of CMC perturbations of the tubes 
TpiT) for all radii p E Ik, when k is large. We shall denote the perturbation 
functions as Wp and $p to emphasize their dependence on p. We shall now 
revert to thinking of these as depending on xq rather than s, and in particular 
we write 

Wp{xo, 6) := wo{xo/p) + w{xo/p, 9). 

Following through the proof, it is not hard to see that these functions depend 
smoothly on p. Furthermore, since the tubes ^(F) already foliate, it suffices to 
verify that the mapping 

(p, xo, 6) ^ G{xo, p(l + wp)e + ^p) (4.33) 

is a local diffeomorphism. 
First, 

\\wp\\L^{SNr) + ||'J'||L~(r,Arr) < cp^. 
Also, from the construction itself, we have 

\\dpWp\\Loc.(^SNr) + \\dp^\\L^{r,Nr) < cp. 

These certainly imply that (|4.33() is a local diffeomorphism, and hence our CMC 
surfaces form a local foliation. □ 



5 Explaining the gaps 

In special cases, such as when F is a circle in the flat torus T"+i — x R"/aZ", 
all of the geodesic tubes about F, of any radius, have constant mean curvature, 
and thus there are no gaps in the local foliation. On the other hand. Theorem 1 
only provides for a local foliation with gaps. As indicated in the introduction, 
there are good reasons why this construction fails from working at all radii. 
We explain this in greater detail now. We first show that for generic metrics, 
the moduli space of CMC surfaces isotopic to a geodesic tubes Tp{T) in M \ F 
is smooth and one-dimensional. The index of the Jacobi operator is constant 
along components, and by estimating the index for the surfaces close to 7^(F), 
p € Ik, we show that there are infinitely many components of this moduli space. 
We conclude by examining in more detail a very degenerate case, where all of 
the geodesic tubes 7^(F) are CMC (when smooth), so there are no gaps, but we 
prove that in this situation there are infinitely many bifurcating branches. 
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5.1 The moduli space 



Denote by Ti.{M, T, g) the moduli space of all CMC surfaces S M which are 
isotopic to any one of the geodesic tubes Tp{T), for p small, in M \ F, with 
respect to the metric g. 

Proposition 5.1 There is an open dense set U of metrics (in the C"^'" topol- 
ogy for any m > 3) on M such that for g Cz U, Ti.{M,T, g) is a smooth one- 
dimensional manifold. 

Fix a surface Eq hi the correct isotopy class, which has CMC with respect to 
some metric go. Nearby surfaces may be written as normal graphs over Sq, and 
hence are parametrized by (small) scalar functions on So- Now consider the 
mapping 

G : C™'"(M, S^T*M) x C^'^CSo) — > ^""^^"(So) 

which assigns to a metric g and a scalar function w on So the mean curvature 
function of the submanifold T,w = {x + w{x)iy{x) : x G Sq}, regarded as a func- 
tion on So, with respect to g. It is not hard to show [Tlj (and also [5]) that the 
differential of this mapping is always surjective, and moreover the restriction of 
this differential to the tangent space of the second factor (which is simply the 
Jacobi operator) is Fredholm of index zero, and is an isomorphism except when 
there exist nontrivial Jacobi fields. The result is then a straightforward appli- 
cation of the Sard-Smale theorem, since G is transverse to the one-dimensional 
curve of constant functions in the range space. 

On the other hand, applying the construction of Theorem 1 when the metric 
g ^ U, we obtain a set of smooth one-dimensional families CMC surfaces in 
T-l{M,T, g). Let Gg{w) = G{g,w)\ an implication of this proposition is that 
when Gg is regular at w and the mean curvature of S^j is equal to TJ, then 
some interval {H — e, H -\- e) parameterizes Ti.{M,T, g) locally near or in 
other words, the mean curvatures of the CMC surfaces near to assume all 
values near to H. Although in this case the corresponding surfaces form a local 
foliation, this need not be true in general. 

5.2 The index 

We now claim that for generic g, 'H{M,r,g) has infinitely many components. 
In the following, for p G Ik, let Ep denote the CMC hypersurface constructed 
in Theorem 1. 

Proposition 5.2 Let g be a metric for which TL{M,T,g) is a smooth one- 
dimensional manifold. Then for each sufficiently large value of k, the surfaces 
Tip, p G Ik, lie in different components of'H{M,T,g) and have index equal to 
Index{T) + 2k + l. 

We prove this theorem by computing the index of T,p\ by definition this is 
the number of negative eigenvalues of the Jacobi operator Cp on Ep. This index 
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is locally constant in H{M,T,g) when this moduh space is nondegenerate, but 
since we shall show that the index increases with k, this will imply the result. 

It follows from H2.26|l and the properties of solution {wq, w, $) given in Propo- 
sition ini that the Jacobi operator about Ep, i.e. the linearization of the operator 
{wq, w, $) I— > npH , has the form 

(5.34) 

where the (linear operator) L satisfies the usual assumptions. 
Lemma 5.1 The quadratic form associated to Cp has the expansion 
Bp{v,^) J^^^{\dsv\^ + \\/gv\^~{n~l)v^) 

+ UJnP /r(|at*|2-(i?(*,Xo)*,Xo)p) 

+ p^C{v,^) 

where v £ UoH^p{SNr) (BUH^piSNT) (i.e. the Sobolev space H^iSNT) defined 
with respect to the vector fields dg and dy.), ^ £ H^{T,NT), a;„ is a positive 
constant depending only on the dimension, and where C is a quadratic form 
satisfying 

\C{v,^)\<c( [ {\dsv\^ + \Vgv\^+v^)+ [ {\dtm^ + \^\^) 
\JSNr Jr 

for some constant c > independent of p. 

Proof : For v and ^' = X^V'j^i i^- this statement, the fibrewise linear 
function on SNT corresponding to 4" is w = We can either regard Bp 

as a quadratic form in the variables (v, ^) or in u + v. Then, by definition 

Bp{v + v)= [ (£p(i;, *))(« + «). 

JSNF 

Inserting the expression H5.34|l and integrating over Sp, we obtain the first 
summand in the expression for Bp, involving only v, without difficulty. Next, 
integrating over the spherical fibres of SNT, we have 

J SNT ^f,^^ J SNT Jr 

there is a similar reduction for the term in ^ of order to an integral over T. 
The error term leading to C arises from the error term in Cp, as well as the 
discrepancy in this last calculation caused by using the volume form on SNT 
rather than the one on Tip. The first of these error terms, 

J L(w,*) (v + v), 
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is almost of the correct form. However, since Cp involves 5"", this error might 
include terms of the form v ^P", which are at first glance too big since, integrating 
by parts, they equal dgV ^' . However, examining the computations leading 
up to (|2.2t)|) . one can check that $" only enters through the term 

p {N, + 0{p^) + L{w, $) + Q{w, $)) 

(recall that we compute p times the mean curvature), which has linearization 

p {N, *'%(! + + L{w, <&) + Q{w, $)) = pi'f", B)p + 0{p^) L{^). 

Hence this gives, at worst, terms like p^ dgV ^' . It is much more straightforward 
to check that all the other terms in C satisfy the correct bounds. □ 

As usual, write u(s, y) — vq{s) + v{s, y), where both sunimands are orthogo- 
nal to the linear eigenfunctions ^j, 1 < j < n; we also identify v(t, y) :— 6)p. 
Thus for V = vo + v + v € UoH^iSNT) QTIH^SNT) © flH^iSNT) we have 
the quadratic form 

Bp{v) =Bo{vo)+p B{v) + B{i) + p2 c{v), 

where 

Bo{vq) := w„ Jpda^uoP - (ri- l)wo) 
B{v) := c^„ (i?(*,Xo)*,Xo)p) 

m ■■= !sr,A\dsi\' + \yev\'-{n-l)v^) 

Assuming that p ^ Ij. = [p'^.p'D, let us now compute the index of B. Since 
this index is locally constant in p, we shall choose 

^/n^ , /I 1 \ 

which is directly in the middle of Ik. 
Clearly, 

j + |Vwp < ci?(w). 

Next, decompose i) — A-v~ , where lies in the sum of the eigenspaces of 
0; with positive or negative eigenvalues, respectively. (Recall that this operator 
is assumed to be nondegenerate, hence has no zero eigenspace.) We then have 

j \dtv+\'' + \v+\^ <cB{v+) 

and 

j + < -cB{v-). 
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Wc can similarly decompose the remaining component f o as Vq +Vq , where 
lie in the eigenspaces corresponding to the positive or negative eigenvalues of 
('Csjvr)o = — {n—1). Using that p lies in the middle of Ik, we can estimate 

I \dsV+\' + \v+\'<cp-'Boiv+) 

and 

J \d,Vo\^ + Ivof < -cp-^SoK")- 
Using all of these estimates, we now obtain that 

\C{v)\ < cp-' (Bo«)-i?o(V)) 

+ c (b{v+) - B{v-)j +cB{v) 

for some constant c > independent of p. This gives, in turn, 

B'{v) < Bp{v) < B"{v) 

where 

B"{v) = {1 + cp) (^Boiv+) + pB{v+)) 

+ {l-cp) (^B{v^) + pB{v-)) 
+ (l + cp2)S({;) 

and 

B'{v) = a-cp) (^Boiv+)+pBiv+)) 

+ {l+cp) (^B{vo)+pB{v-)) 
+ {l-cp'')B{v). 
These upper and lower bounds on B imply that 

Index (B") < Index (B^) < Index (B'). 
The proof is completed by the following 

Lemma 5.2 When p is small enough, the index of B' and B" are both equal 
to Index{T) + 2k + l. 

Proof : If /9 is chosen so that l — cp> 1/2, then the index of B' and B" equals 

the sum of the dimensions of the spaces on which Bq and B are negative. But 
these equal 2k + 1 and Index (F), respectively. □ 
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5.3 Bifurcations in a degenerate case 

We consider in more detail the (very) degenerate case where (M, g) is the flat 
torus T"+i = 51 X ]R"/aZ" for a sufficiently large, and T = x {0}. After 
moding out by all the continuous symmetries, the moduli space is still one- 
dimensional, but has infinitely many singularities. 

Each of the geodesic tubes 7;,(rj := = 5^ x ^"-^(p) is CMC, with mean 
curvature Hp = {n — l)/p. We only consider the case where Sp is embedded, 
i.e. when p < 1/2. The Jacobi operator for is 

Lp = Aj:^ + \Aj:f. 

In terms of our standard cylindrical coordinates, t G S^, 6 G 

Introducing eigendata {(j)k{t),k^} and {ipi{9),—\'j) in each component, we see 
that Lp reduces to multiplication by B{k, £, p) = —k"^ + p~'^{n — 1 — A^) on the 
(fc, £) eigenspace. Since Xj is always of the form — 2+j) for some nonnegative 
integer j, and hence n — 1 — < unless £ — 0. This gives the 

Proposition 5.3 The surface Sp is always degenerate; its nuUspace consists of 
the span of the eigenmodes (j)o{t)'ipi{9), i= 1, . . . ,n and, in case = (n— 1)/A;^ 
for some A: S N, also (f>k{t)tpa{ff}- 

The 'trivial degeneracies' are those comprised by the first set of elements, which 
exist for all p. These correspond to the obvious geometric fact that translating 
Ep parallel to itself in any direction normal to F gives a family of CMC surfaces 
with the same mean curvature. These can be eliminated if we mod out by these 
symmetries. Namely, using linear coordinates x = {xo,xi, . . . ,a;„) in T"+^, let 
G be the finite group generated by reflections in the xj = plane, j = 1, . . . ,n. 

Coroll£iry 5.1 Acting on the space of G -invariant functions on SNT, the Ja- 
cobi operator Lp is degenerate if and only if p^ = {n — l)/k^, k G N. 

These degeneracies have a direct geometric explanation too, for it is precisely at 

these radii where the family of 'fc-bump' Delaunay surfaces begins to develop. 

We now have the picture that H{T^~^^ , , go) consists of the union of the 
interval (0, a/2)p and infinitely many other intervals (0, {n—l)/k'^]e,,, where the 
variable et is the Delaunay necksize in the k^^ bifurcating branch. In other 
words, the moduli space looks like an open interval with infinitely many spines 
sticking out of it. When the metric on T"+^ is perturbed generically, this moduli 
space smooths out; the singularity at each degenerate radius disappears, and 
this 'spiny interval' breaks into infinitely many components. The CMC surfaces 
for these slightly perturbed metrics are small perturbations either of the geodesic 
tubes or else of the Delaunay surfaces, except near the turning points. The gaps 
encountered in our construction correspond exactly to the small regions around 
these degenerate radii where the moduli space is curving away from the interval 
(0,a/2)p. 
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6 Limits of constant mean curvature foliations 



There is a sort of converse to Theorem 11.11 which we can prove regardless of 
the dimension of F. Let F be a closed ^-dimensional submanifold of M"^^, 
1 < £ < n, and suppose that there exists a sequence of hypersurfaces Sj such 
that 

a) has constant mean curvature (n — £)/npj, where {pj} is a decreasing 
sequence with pj \0; 

b) is isotopic in M \ F to the tube Tp^ (F) ; 

c) There exists a c > such that Ej is contained inside Tcp- (F) for aU j. 

Item c) implies that Ej F in Hausdorff distance; note also that we are not 
requiring the existence of any sort of local foliation, just a sequence of CMC 
hypersurfaces converging to F. We conjecture that some set of hypotheses very 
near to these (for example, assuming also a bound on the second fundamental 
form, as in the introduction, or that E-,- is trapped between the tubes of radius 
c' Pj and cpj for fixed constants < c' < c.) should be enough to ensure that 
F is minimal. Unfortunately, this seems to be quite difficult to prove, and so 
we shall restrict ourselves to a very special situation by making a fourth, quite 
restrictive, hypothesis: 

d) There exist functions Wj G €^'"{3 NT), $j e C'^'"{r,NT) satisfying 

\\Wj\\c2,c + ||$j||c2,o < 

for some c > 0, independent of pj, such that 

Ej = Tp^{wj,<i>j). 

Quite important (and potentially restrictive) here is that the norms of Wj and 
$j are bounded in C^'", not Cp'". 

Theorem 6.1 Let F be a compact embedded (.-dimensional submanifold of 
M for some 1 < £ < n, and suppose that Ej is a sequence of CMC hypersurfaces 
converging to F and satisfying the hypotheses a) - d). Then F is minimal. 

The proof is based on an argument from geometric measure theory which is 
now fairly standard in the analysis of such 'condensation problems', cf. [2]. 

We drop the subscript j and consider a functional for which each Ep is 
critical. The argument proceeds by writing the formula which expresses the 
fact that the first variation of this functional vanishes, and then taking the limit 
of this formula (in a very weak sense) as p ^ 0. The limiting first variation 
equation implies the minimality of F. 

Any one of the CMC hypersurfaces E = Ep bound a compact domain £'(E), 
which is the component of M \ E containing F. Define the measure 

d/is = dA^ — nH dV^: , 
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where H is the (constant) mean curvature of S and where we have set 

dA^ rr Ls and dV^ := 

(Ti*^ is k dimensional Hausdorff measure). 

CMC hypersurfaces are critical for the functional S / dji^. In other 
words, if X is any vector field on A/, and the associated one-parameter 
family of diffeomorphisms, then 



d^l^lT. 



= 0. 



(6.35) 



We now compute this variation another way. In fact, for any hypersurface E 
and any continuous function /, we derive that 

dt J fdA^.^ = J XfdA^+ J /(divX - (VA,X,iV))dAs, 

where N is the unit normal to S, and similarly. 



dt J fdV^.^ 



= / XfdVj: 



fdivXdVs. 



Hence, setting / = 1 and using H6.35|l . we obtain 

= y /d^i^js = J diyX dfi^ - y (VjvX, N) dA^. 



(6.36) 



Next, let dLr = Ti.^ [r and denote by Un-e the volume of S" ^ with respect 
to its standard metric (thus LUn-i/{n + 1 — ^) is the volume of the ball B"''"^^^). 



Lemma 6.1 j4s p \ 0, 

w„_^ fiLr, and p"- 
in the sense of measures. In particular, 



P dpy. — — : dL-i 



n + l-i 



?i + 1 - . 



■dLr 



Proof : This follows from Fubini's theorem and the fact that the functions w 
and $ appearing in the parametrization of E are uniformly controlled in C^'°' 
as p ^ 0. □ 

On the other hand, we also have 
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Lemma 6.2 Let Ei, . . . , En+i-i he a local orthonormal frame for NT. Then 
as p\0, 

p^-" {VNX,N}pdAs^ - y^{\/E^X,E,)pdLr 

n + 1 — I 

in the sense of measures. 

Proof : As before, using Fubini's theorem and the uniform control on w and 
$, it suffices to check that this formula holds for the round sphere of radius p 
in K"+^ as p — + 0, where again the formula is a standard computation. □ 

Now multiply Ht).36|l by p^"" and let p ^ 0. From these two lemmas, we 
conclude that 

dhfX ^^{VE^X.Ei)^ dLr = 

On the other hand, if i^i, . . . , i^^ is a local orthonormal frame for TF, then 

i n+i-e 
divX = ^(VF^X,F,)+ ^ {^E,X,E,), 
j=i i=i 

so this last equation is equivalent to 

J Y,{X,\/F^F,)dLr^O. 

3 

Since the vector field X is arbitrary, we conclude that the normal component 
of \I FjFj is equal to 0. This implies that the mean curvature of F vanishes, 
i.e. that F is minimal. □ 

As already discussed at the beginning of this section, it would be much 
nicer to prove this theorem under less stringent hypotheses. We conclude by 
discussing this in more detail. 

Suppose that there exist sequences of intervals Ik — {p'ki P'k) with, with 

P'k^P'k ~^ 0' such, that for each k and p € Ik there exists a CMC hypersurface Ep 
isotopic to ^(F) in M\F. Suppose furthermore that this hypersurface satisfies: 

a') the hypersurfaces {Epjpg/j, form a local foliation; 

b') the mean curvature of Ep equals 

c') there exists a constant c > 0, independent of k and p g /fc, such that 
Ep C Bp{T) := {q e A/"+i : dist<,((7,F) < cp}. 

d') There exists a constant c > 0, again independent of k and p & Ik, such 
that < c/p. 
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We conjecture that these hypotheses alone are sufficient to conclude that F is 
minimal. Indeed, it is possible to prove many of the necessary facts, but a few 
crucial ones seem much more difficult to obtain. 

If we rescale Ep from a point p g F by the factor 1/p, we obtain a family 
of CMC surfaces which are cylindrically bounded. We can obtain area bounds 
for these rescaled hypersurfaces, just as in m, and so conclude that at least 
along subsequences, this family converges to a complete embedded cylindrically 
bounded hypersurface in It is known (and also that all such 

hypersurfaces must lie in the family of Delaunay unduloids D^; however, using 
that the Ep are leaves of a local CMC foliation, we obtain a global bounded 
positive Jacobi field on the limiting surface, and it may be checked directly that 
in the Delaunay family, only the cylinder admits such a Jacobi field. Thus far 
we have shown that, having fixed p G T and rescaling about this point, then 
some subsequence of these surfaces, say Ej, converges to a cylinder with axis 
parallel to the rescaled limit of F. The first difficulty is a standard one in the 
subject: if we knew the uniqueness of this limit, then we could straightaway 
conclude the existence of functions wj and such that Ej = Tp- {wj,^j). We 
would also be able to conclude that 



on the other hand, it is quite straightforward to prove that Jjjf^^ ^ dVs < 
(jp^+^-t _ In any event, we are only able to show that the second conclusion of 
Lemma 16.11 holds, and so we are unable to take the weak limit of d/is^ . The 
final difficulty arises in proving the analogue of Lemma [6.21 and this is the case 
because of the rather weak control we have for the derivatives of w and $ in the 
F direction. 

Plausibly, the most general theorem of this sort would involve a sequence 
of CMC hypersurfaces which are known to converge in Hausdorff distance, and 
perhaps satisfying hypotheses a') - d'). The conclusion should be that Ep 
converge to a minimal submanifold F, of some dimension £, away from a set of 
Hausdorff measure smaller than £. For example, if £ = 1, it is quite conceivable 
that such a set of surfaces might converge to a broken geodesic (satisfying certain 
constraints at the break points). 
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